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The relativistic theory of hydrodynamic fluctuations, or noise, is de- 
rived and applied to high energy heavy ion collisions. These fluctuations 
are inherent in any space-time varying system and are in addition to ini- 
tial state fluctuations. We illustrate the effects with the boost-invariant 
Bjorken solution to the hydrodynamic equations. Long range correlations 
in rapidity are induced by propagation of sound modes. The magnitude of 
these correlations are directly proportional to the viscosities. These fluctu- 
ations should be enhanced near a phase transition or rapid crossover. 

PACS numbers: 25.75.Ag, 25.75.Gz, 25.75.Ld 

1. Introduction 

Cosmology has experienced tremendous advances in the past decade. 
These advances have been driven by observations. Much of the information 
comes from fluctuations in temperature of the cosmic microwave background 
radiation (CMBR) as observed by the Wilkinson Microwave Anisotropy 
Probe (WMAP) [I]. The results are usually expressed in terms of the power 
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spectrum, 1(1 + l)Ci/2ir versus the multipole moment I, where 

T{0,<l>) = $> Im Y lm (M) 
Im 

Q = (\a lm \ 2 ) (1) 

and where the averaging is done over points in the sky. In WMAP7 the 
power spectrum extends up to values of I on the order of 1000. At the 
Relativistic Heavy Ion Collider (RHIC) and at the Large Hadron Collider 
(LHC) jets clearly stand out against a background of low transverse mo- 
mentum particles in a lego plot. In Pb-Pb collisions at the LHC these jets 
can have transverse momentum as high as 100 to 150 GeV/c compared to 
the background particles of several GeV/c [2]. Here we are not interested in 
the jets but in the fluctuations in the background. What information can 
we extract from these fluctuations? 

The theory behind fluctuations in the CMBR is highly nontrivial |3] 
even if the basic ideas are rather intuitive. One should expect a similar 
nontriviality in the theory for high energy heavy ion collisions. There are 
at least four important sources of fluctuations in such collisions: 

(i) Initial state fluctuations: These arise because nuclei are composed 
of nucleons which in turn are composed of quarks and gluons. The 
fluctuations are both statistical and quantum in nature. 

(ii) Hydrodynamic fluctuations: These arise due to finite particle number 
fluctuations in a given coarse-grained fluid cell. 

(iii) Fluctuations induced by jets: Jet production is a random process. 
They deposit energy and momentum as they traverse the system. 

(iv) Freeze- out fluctuations: These arise when converting a coarse-grained 
fluid cell into individual particles which are subsequently described by 
a Monte Carlo transport model. 

The goal here is to develop a relativistic theory of hydrodynamic fluctua- 
tions for application to high energy heavy collisions. After describing the 
basic theory we will apply it to the boost-invariant solution of Bjorken for 
high energy collisions. Although not realistic enough to compare directly 
with data it will demonstrate that correlations develop whose magnitude 
and shape are controlled, in a quantitative manner, by the shear and bulk 
viscosities as well as by the equation of state. 

Intuitively, hydrodynamic fluctuations or noise becomes important when 
gradients of temperature, density, and particle composition become large. 
In such situations the coarse-grained fluid cells must be relatively small to 
adequately represent these gradients. The smaller the cell the more impor- 
tant fluctuations will be. Certainly this is the case with high energy heavy 
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ion collisions. Analogous situations may be found throughout the physical, 
chemical and biological literature; see [1] and citations within. For example, 
Eggers [5 J performed a theoretical study of the breakup of liquid nanojets 
with the conclusion that "noise is the driving force behind pinching, speed- 
ing up the breakup to make surface tension irrelevant" . Similar conclusions 
were reached by Kang and Landman [6] who studied the breakup of liquid 
nanobridges with a molecular dynamics approach, with a lubrication equa- 
tion (smooth fluid dynamics), and with a stochastic lubrication equation. 
Inclusion of noise in the lubrication equation provided results very similar 
to the molecular dynamics simulations. 

With these as motivations we proceed to the general study of relativistic 
hydrodynamic fluctuations in the next section followed by an application to 
the boost invariant hydrydynamics of Bjorken. 

2. Relativistic Hydrodynamic Fluctuations 

Now we turn to the topic of hydrodynamic fluctuations. The energy- 
momentum tensor density for a perfect fluid is 

2£L = + Wu" • (2) 

Here w = P + e = Ts + \in is the local enthalpy density, fi is the baryon 
chemical potential, n is the baryon density, and u M is the local flow velocity. 
The metric is (+, —,—,—). Corrections to this expression are proportional 
to first derivatives of the local quantities whose coefficients are the shear 
viscosity 77, bulk viscosity £, and thermal conductivity x- Explicit expres- 
sions may be found in textbooks [8] which are useful to summarize here. 
Dissipative contributions are added to the energy-momentum tensor and 
baryon current as follows. 

J* = nu^ + AJV (3) 

There are two common definitions of the flow velocity in relativistic dissi- 
pative fluid dynamics which are important to distinguish. In the Landau- 
Lifshitz approach is the velocity of energy transport. In the Eckart 
approach u M is the velocity of baryon number flow. In high energy heavy 
ion collisions, at the upper range of RHIC energies and at the LHC, the 
net baryon number is very small compared to the entropy density or to 
the number of baryons plus anti-baryons. Therefore the Landau-Lifshitz 
approach is the relevant one. 

In the Landau-Lifshitz approach the dissipative part of the energy- 
momentum tensor satisfies u (i AT' li/ = 0. The most general form of the 
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energy-momentum tensor is 

AT^ = v (A»u u + A V) + (§77 - C) H^d -u. (4) 

Here 

= u»u u - (5) 

is a projection tensor normal to 

A M = df, - Uf, (u ■ d) (6) 

is a derivative normal to u^, and 

Q a = d a T-T{u- 0) u a (7) 

is a heat flow vector whose nonrelativistic limit is Q = — VT. The baryon 
current is modified by 

AJM = x (v) 2 AM(/3/i) ' (8) 

which satisfies u^AJ^ = 0. This insures that n is the baryon density in the 
local rest frame. The entropy current in this approach is 

s m = su » - (3fiAJ^ . (9) 
In the local rest frame entropy is generated according to the divergence 



+ |(V-u) 2 + ^(VT + Tu) 2 . (10) 

The term Tu is a relativistic correction to VT, being smaller by a factor 
of 1/c 2 in physical units. All three dissipation coefficients must be non- 
negative to insure that entropy can never decrease. 

Next we add small fluctuations to the energy-momentum tensor 

= T^ + AT< XV + (11) 

and to the baryon current 

= n / + AJ» + J" . (12) 

The fluctuations must satisfy the conditions u^S^ = and u^I^ = 0, 
which are the same conditions satisfied by AT^ and by AJ. These fluctu- 
ations have zero average value at every space-time point. The averaging is 
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done with an ensemble of heavy ion collisions, all prepared with exactly the 
same initial conditions when the hydrodynamic description may be applied. 
However, the average of a product of fluctuations is not necessarily zero. 
We follow section 88 on hydrodynamic fluctuations of [9] to derive these 
correlators. After some analysis we find 



(S^(xi)S Q/3 (x 2 )> = 2T [rj (H' ux H uf> + H^H 1 

+ (C - §77) H^H al3 ] 6{xi - x 2 ) (13) 



and 



{I»{x x )I v {x 2 )) = 2x (~) H^5( Xl -x 2 ). (14) 

The mixture {S^ '{x\)I a (x 2 )) naturally enough is zero. These correlation 
functions have their origin in the fluctuation-dissipation theorem. There are 
two essential observations concerning them. First, they are proportional to 
Dirac delta-functions. Different coarse-grained fluid cells are assumed to be 
independent. Second, the magnitudes of the correlations are directly pro- 
portional to the shear and bulk viscosities and to the thermal conductivity. 
This is where the microscopic physics lies. 

The procedure for implementing these stochastic sources are as follows. 

• Solve the hydrodynamic equations for an arbitrary source function. One 
can imagine doing this in principle although it may be difficult in practice. 

• Perform averaging using the above expressions for the sources to obtain 
observable correlation functions. 

• The stochastic fluctuations may or may not be perturbative, depending 
on the physical conditions. 

To develop an understanding of the implementation we consider a simple 
model for heavy ion collisions in the next section. 



3. Boost Invariant Model 

It is insightful to work out a particular example; we shall do this for 
the well-known boost-invariant Bjorken model. In the absence of fluctua- 
tions the important results are that the temperature depends only on the 
proper time as T = T(t), the entropy density decreases with proper time 
as s(t) = soTq/t, where tq is the equilibration time, and the fluid velocity 
is = (cosh£, 0, 0, sinh£), where £ is the space-time rapidity. Space-time 
fluctuations may be expressed as 



T = T(t)+5T(£,t) 

u 11 = (cosh(£ + w(£,r)),0,0,sinh(£ + a;(e,r))) 



(15) 
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The function oj is dimensionless, and it is convenient to also use the di- 
mensionless variable p = Ss/s instead of ST (the latter are related by ther- 
modynamic identities). For purposes of illustration, these fluctuations are 
treated as perturbations. After lengthy calculation one finds the typical 
linear response relations 

f T dr' ~ 

X(k,r) = - G x (k;r,T')f(k,T') (16) 

Jto T 

where X is either p ox uj. This relation is given in terms of the variable k 
which corresponds to the Fourier transform of the variable £. The function 
/ is the single scalar function representing the noise which can expressed as 



The correlation functions are 

r T f dr T(t) 



(17) 



w^/)no,T,)> = ^jr 



T 3 W 2 (t) 



Hr) + C(r) 



/oo 
dke lk ^G X Y(k;Tf,T) . (18) 
-oo 

Here GxY(k;Tf,T) = r)GV(— £;; 7y , r), A is the effective trans- 

verse area of the colliding nuclei, and tj is the freeze-out time at which the 
transition from hydrodynamic flow to free-streaming of particles takes place. 
The Green functions Gx(k,T,T') are even functions of k and calculable. 

Take, for example, the equation of state P = |e, and treat the underlying 
expansion as essentially inviscid. Then 



GJk;r,r')= - 



2 + 37 - 97 2 / r y 2 - 37 - 97 2 / r 



67 \t' J 67 Vt', 

(19) 

where 7 = — 3k 2 . Note that 7 may be real or imaginary. Similar ex- 
pressions can be written down for the other response functions. The Fourier 
transformed functions are singular, with Dirac delta-functions and deriva- 
tives of them at £ = and at the sound horizon £ = 2v s In (r/r'). Figure 1 
shows the regular part of G pp while Fig. 2 shows the singular part which is 
smeared by a Gaussian of arbitrary with to display the singularities. The 
origin of the singularities are the space-time delta-functions in the original 
correlation functions; this represents white noise. It is possible to cure these 
singularities by using finite range correlations, which is colored noise. 
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Fig. 1. The regular (continuous) part of the correlator G pp (£;t/,t) with v 2 = 1/3 
and ln(r//r) = 4. Note the sound horizon at £ = 2v s ln(r//r). 




Fig. 2. The singular part of the correlator G pp (£;t/,t) with v 2 = 1/3 and 
ln(r//r) = 4. The function is smeared by a Gaussian of variance a 2 — 0.1 in 
order to show the nature of the singularities. 
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Fluctuations in the local temperature and flow velocity fields give rise 
to a nontrivial 2-particle correlation when the fluid elements freeze-out to 
free-streaming hadrons. In momentum-space rapidity this correlation is 



/ dNfa) dN(rn) _ /^\ 2 \ /^lY 1 = 

\ drj dr\ \ dr\ / / \ drj / 7r 4 A r e fr TfTf 




I K(A V ) 
o 

(20) 

where N e s is the effective number of bosonic degrees of freedom at the 
initial time To and temperature To and d s is the spin/isospin degeneracy of 
the hadron species. The correlation is directly proportional to the ratio of 
shear viscosity to entropy density, (77/5)0, assumed here to be temperature 
independent. The function K(Ar]) can be computed numerically for a given 
hadron mass. It is shown in Fig. 3 for the choice To = 600 MeV, Tj = 150 
MeV, tj = 10 fm, and N e g = 47.5. Folding the fluctuations with the 
thermal distribution function smooths out the singularities. Incorporation 
of viscosity in the expansion dynamics and finite range correlations fill-in the 
dip at A77 = 1.5 but otherwise does not much affect the shape or magnitude 
of K. It reminds one very much of the near-side ridge |10[ [Ti] \V2\. 




Fig. 3. The correlation function K (A77) in the dN/drj fluctuations. 
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4. Conclusion 

In this work we have derived the general theory of relativistic hydrody- 
namic fluctuations and showed how to apply them to high energy heavy ion 
collisions. We illustrated the general procedure with the simple example 
of boost-invariant hydrodynamics. This example suggests long-range cor- 
relations in rapidity caused by hydrodynamic fluctuations in temperature 
and flow velocity. Allowing for transverse expansion should also give rise to 
correlations in azimuth (13] . The magnitude of the correlations are directly 
proportional to the viscosities, and the range and shape are sensitive to the 
speed of sound of the medium. Many more details and references are to be 
found in an upcoming publication. 

To fully explore the implications of the ideas and formulas presented 
here and for detailed comparison to data requires a full 3+1 dimensional 
relativistic dissipative hydrodynamic code such as in [II]. It will also be 
interesting to explore the implications for a critical point in the QCD phase 
diagrams where fluctuations should be enhanced [15J. 

Our conclusion is that fluctuations are interesting and can provide im- 
portant information on transport coefficients. We are learning, we are en- 
joying, and there is plenty of work ahead for both theorists and experimen- 
talists. 

Finally I congratulate Johann Rafelski on his 60th birthday! I have 
known Jan for more than 30 years, but we have collaborated on only one 
project. With Berndt Muller we put together an annotated reprint collec- 
tion entitled Quark- Gluon Plasma: Theoretical Foundations. We assembled 
what we considered to be the pioneering papers in this field in 13 different 
categories. The oldest paper reprinted was Fermi's 1950 article. The newest 
papers were from 1992-93 after which essentially all articles became easily 
accessible and searchable on the preprint archive. The collection is over 800 
pages, and was published by Elsevier in 2003. I highly recommend it for 
newcomers for a primer on the history of our field. Thank you, Jan, for 
having the original idea, and for including me in that project! 
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